Overview
This document is meant to serve as an update/check-in for the “Snow Fire” project, as well as just a way for me to organize my thoughts. Previously I was approaching this project using the more available 50m resolution snow water equivalent data collected by ASO. However, I’m switching to a more fine-scale approach following feedback from others and some concerns about biased results using the coarser data.
At this point, this new approach is heavy on the thinking out loud and thought-organization, but in the future (and once we get all the data), I’ll focus more summarizing analysis and results. Please let me know if you have any comments/suggests.
Motivation
Snow depth varies across forested areas, in part as a function of tree cover, which influences both snow interception and ablation. This likely varies both under canopy vs. in forest gaps, but also with gap internal distance from gap edge, and consequently with gap size and shape.
Broad management strategies have resulted in differences in forest coverage and configuration across the region. Specifically, timber harvest history and fire suppression vs. managed wildfire, determine the size, shape and location of forest gaps (as well as the size and degree of shading of edge trees).
Fire regimes are changing with burns creating increasingly large forest gaps. Also fires at higher elevations are becoming more common.
In contrast to reference conditions, within fire-suppressed landscapes gaps are likely fewer and smaller where fire has been absent, and larger where recent fires have occurred.
The climate is changing which affects both the magnitude of the snow pack at different elevations and timing of snow-melt. Also climate change is leading to higher densities of high-elevation forest.
The amount of snow and the timing of its melt is important for a lot of things like life, almonds, and swimming pools.
Data Availability
We are using Nasa’s 3m lidar-derived ASO data. Lidar fights have occurred starting in 2013, although they improved their equipment in 2014, so let’s start that year. They have conducted flights across a number of basins, but we’ll focus on the Merced and Kings basins since these have had the closest to natural fire regimes and likely most ‘natural’ distributions of gap types. If we want more temporal and/or spatial replication we could look at other basins, but we would need to wrangle a lot of data, and data requests haven’t been rapid so far.
Currently we have data for 2014 from the Merced basin, and a request in for all other flights in the Merced and Kings basins.
d <- read.csv("data/reference/snowdepth3m_avail.csv") %>%
mutate(date = as.character(date),
date = as.Date(date, "%Y%m%d"),
year = year(date), #year
jday = as.integer(format(date, '%j')))
## Plot flights across years and day
ggplot(d, aes(x = jday, y = year, color = basin, shape = factor(pocess))) +
geom_point(size = 4) +
scale_shape_manual(values = c(21, 19),
labels = c("Requested","Obtained"),
name = "Data Status") +
scale_color_brewer(palette = "Dark2",
labels = c("Kings", "Merced"),
name = "Basin") +
xlab("Julian Day") + ylab(NULL)

Questions
There are (at least) two ways to approach this analysis using 3m ASO data:
At the gap level: How does gap size and shape influence mean snow depth and snow disappearance date? Could apply the SDC metric here as well.
At the pixel level: How does snow depth and presence vary with distance to patch edge? Do pixels in shaded edges hold snow longer?
Either way we’d want to incorporate elevation, aspect, date, and snow-year to ask the follow if possible:
- Do this relationship depend upon climate (both inter-annual and across elevation)?
Relevant Literature
There are clearly more than is listed below. If there are any important omissions at this point, please let me know.
Related to fire effects on forest gaps/patches: (Steel, Koontz, and Safford 2018; Stevens et al. 2017; Collins et al. 2017; Boisram’e et al. 2017)
Lidar measurements of snow depth: (A. A. Harpold et al. 2014; Painter et al. 2016)
Snow and tree cover: (Adrian A. Harpold et al. 2014; Currier and Lundquist 2018; Lundquist et al. 2013; Stevens 2017)
Bibliography
Boisram’e, Gabrielle, Sally Thompson, Brandon Collins, and Scott Stephens. 2017. “Managed Wildfire Effects on Forest Resilience and Water in the Sierra Nevada.” Ecosystems 20 (4): 717–32. doi:10.1007/s10021-016-0048-1.
Collins, Brandon M., Jens T. Stevens, Jay D. Miller, Scott L. Stephens, Peter M. Brown, and Malcolm P. North. 2017. “Alternative Characterization of Forest Fire Regimes: Incorporating Spatial Patterns.” Landscape Ecology 32 (8): 1543–52. doi:10.1007/s10980-017-0528-5.
Currier, William Ryan, and Jessica D. Lundquist. 2018. “Snow Depth Variability at the Forest Edge in Multiple Climates in the Western United States.” Water Resources Research 54 (11): 8756–73. doi:10.1029/2018WR022553.
Harpold, A. A., Q. Guo, N. Molotch, P. D. Brooks, R. Bales, J. C. Fernandez-Diaz, K. N. Musselman, et al. 2014. “LiDAR-Derived Snowpack Data Sets from Mixed Conifer Forests Across the Western United States.” Water Resources Research 50 (3): 2749–55. doi:10.1002/2013WR013935.
Harpold, Adrian A., Joel A. Biederman, Katherine Condon, Manuel Merino, Yoganand Korgaonkar, Tongchao Nan, Lindsey L. Sloat, Morgan Ross, and Paul D. Brooks. 2014. “Changes in Snow Accumulation and Ablation Following the Las Conchas Forest Fire, New Mexico, USA.” Ecohydrology 7 (2): 440–52. doi:10.1002/eco.1363.
Lundquist, Jessica D., Susan E. Dickerson-Lange, James A. Lutz, and Nicoleta C. Cristea. 2013. “Lower Forest Density Enhances Snow Retention in Regions with Warmer Winters: A Global Framework Developed from Plot-Scale Observations and Modeling.” Water Resources Research 49 (10): 6356–70. doi:10.1002/wrcr.20504.
Painter, Thomas H., Daniel F. Berisford, Joseph W. Boardman, Kathryn J. Bormann, Jeffrey S. Deems, Frank Gehrke, Andrew Hedrick, et al. 2016. “The Airborne Snow Observatory: Fusion of Scanning Lidar, Imaging Spectrometer, and Physically-Based Modeling for Mapping Snow Water Equivalent and Snow Albedo.” Remote Sensing of Environment 184 (October): 139–52. doi:10.1016/j.rse.2016.06.018.
Steel, Zachary L., Michael J. Koontz, and Hugh D. Safford. 2018. “The Changing Landscape of Wildfire: Burn Pattern Trends and Implications for California’s Yellow Pine and Mixed Conifer Forests.” Landscape Ecology 33 (7): 1159–76. doi:10.1007/s10980-018-0665-5.
Stevens, Jens T. 2017. “Scale-Dependent Effects of Post-Fire Canopy Cover on Snowpack Depth in Montane Coniferous Forests.” Ecological Applications 27 (6): 1888–1900. doi:10.1002/eap.1575.
Stevens, Jens T., Brandon M. Collins, Jay D. Miller, Malcolm P. North, and Scott L. Stephens. 2017. “Changing Spatial Patterns of Stand-Replacing Fire in California Conifer Forests.” Forest Ecology and Management 406 (December): 28–36. doi:10.1016/j.foreco.2017.08.051.
LS0tCnRpdGxlOiAiU25vd0ZpcmUgVXBkYXRlIDIiCmF1dGhvcjogIlphY2sgU3RlZWwiCmRhdGU6ICdgciBmb3JtYXQoU3lzLkRhdGUoKSwgIiVZLSVtLSVkIilgJwpiaWJsaW9ncmFwaHk6IFsuLi8uLi9yZWZlcmVuY2VzLmJpYl0KIyBvdXRwdXQ6CiMgICBodG1sX25vdGVib29rOgojICAgICBjb2RlX2ZvbGRpbmc6IGhpZGUKIyAgICAgdG9jOiBubwpvdXRwdXQ6IHJtYXJrZG93bjo6Z2l0aHViX2RvY3VtZW50Ci0tLQoKIyBPdmVydmlldyAKClRoaXMgZG9jdW1lbnQgaXMgbWVhbnQgdG8gc2VydmUgYXMgYW4gdXBkYXRlL2NoZWNrLWluIGZvciB0aGUgIlNub3cgRmlyZSIgcHJvamVjdCwgYXMgd2VsbCBhcyBqdXN0IGEgd2F5IGZvciBtZSB0byBvcmdhbml6ZSBteSB0aG91Z2h0cy4gUHJldmlvdXNseSBJIHdhcyBhcHByb2FjaGluZyB0aGlzIHByb2plY3QgdXNpbmcgdGhlIG1vcmUgYXZhaWxhYmxlIDUwbSByZXNvbHV0aW9uIHNub3cgd2F0ZXIgZXF1aXZhbGVudCBkYXRhIGNvbGxlY3RlZCBieSBBU08uIEhvd2V2ZXIsIEknbSBzd2l0Y2hpbmcgdG8gYSBtb3JlIGZpbmUtc2NhbGUgYXBwcm9hY2ggZm9sbG93aW5nIGZlZWRiYWNrIGZyb20gb3RoZXJzIGFuZCBzb21lIGNvbmNlcm5zIGFib3V0IGJpYXNlZCByZXN1bHRzIHVzaW5nIHRoZSBjb2Fyc2VyIGRhdGEuICAKCkF0IHRoaXMgcG9pbnQsIHRoaXMgbmV3IGFwcHJvYWNoIGlzIGhlYXZ5IG9uIHRoZSB0aGlua2luZyBvdXQgbG91ZCBhbmQgdGhvdWdodC1vcmdhbml6YXRpb24sIGJ1dCBpbiB0aGUgZnV0dXJlIChhbmQgb25jZSB3ZSBnZXQgYWxsIHRoZSBkYXRhKSwgSSdsbCBmb2N1cyBtb3JlIHN1bW1hcml6aW5nIGFuYWx5c2lzIGFuZCByZXN1bHRzLiBQbGVhc2UgbGV0IG1lIGtub3cgaWYgeW91IGhhdmUgYW55IGNvbW1lbnRzL3N1Z2dlc3RzLiAgCgojIE1vdGl2YXRpb24gIAoKIDEuIFNub3cgZGVwdGggdmFyaWVzIGFjcm9zcyBmb3Jlc3RlZCBhcmVhcywgaW4gcGFydCBhcyBhIGZ1bmN0aW9uIG9mIHRyZWUgY292ZXIsIHdoaWNoIGluZmx1ZW5jZXMgYm90aCBzbm93IGludGVyY2VwdGlvbiBhbmQgYWJsYXRpb24uIFRoaXMgbGlrZWx5IHZhcmllcyBib3RoIHVuZGVyIGNhbm9weSB2cy4gaW4gZm9yZXN0IGdhcHMsIGJ1dCBhbHNvIHdpdGggZ2FwIGludGVybmFsIGRpc3RhbmNlIGZyb20gZ2FwIGVkZ2UsIGFuZCBjb25zZXF1ZW50bHkgd2l0aCBnYXAgc2l6ZSBhbmQgc2hhcGUuICAKIAogMi4gQnJvYWQgbWFuYWdlbWVudCBzdHJhdGVnaWVzIGhhdmUgcmVzdWx0ZWQgaW4gZGlmZmVyZW5jZXMgaW4gZm9yZXN0IGNvdmVyYWdlIGFuZCBjb25maWd1cmF0aW9uIGFjcm9zcyB0aGUgcmVnaW9uLiBTcGVjaWZpY2FsbHksIHRpbWJlciBoYXJ2ZXN0IGhpc3RvcnkgYW5kIGZpcmUgc3VwcHJlc3Npb24gdnMuIG1hbmFnZWQgd2lsZGZpcmUsIGRldGVybWluZSB0aGUgc2l6ZSwgc2hhcGUgYW5kIGxvY2F0aW9uIG9mIGZvcmVzdCBnYXBzIChhcyB3ZWxsIGFzIHRoZSBzaXplIGFuZCBkZWdyZWUgb2Ygc2hhZGluZyBvZiBlZGdlIHRyZWVzKS4gIAoKICAgIGEpIEZpcmUgcmVnaW1lcyBhcmUgY2hhbmdpbmcgd2l0aCBidXJucyBjcmVhdGluZyBpbmNyZWFzaW5nbHkgbGFyZ2UgZm9yZXN0IGdhcHMuIEFsc28gZmlyZXMgYXQgaGlnaGVyIGVsZXZhdGlvbnMgYXJlIGJlY29taW5nIG1vcmUgY29tbW9uLiAgCiAKICAgIGIpIEluIGNvbnRyYXN0IHRvIHJlZmVyZW5jZSBjb25kaXRpb25zLCB3aXRoaW4gZmlyZS1zdXBwcmVzc2VkIGxhbmRzY2FwZXMgZ2FwcyBhcmUgbGlrZWx5IGZld2VyIGFuZCBzbWFsbGVyIHdoZXJlIGZpcmUgaGFzIGJlZW4gYWJzZW50LCBhbmQgbGFyZ2VyIHdoZXJlIHJlY2VudCBmaXJlcyBoYXZlIG9jY3VycmVkLiAgCgogMy4gVGhlIGNsaW1hdGUgaXMgY2hhbmdpbmcgd2hpY2ggYWZmZWN0cyBib3RoIHRoZSBtYWduaXR1ZGUgb2YgdGhlIHNub3cgcGFjayBhdCBkaWZmZXJlbnQgZWxldmF0aW9ucyBhbmQgdGltaW5nIG9mIHNub3ctbWVsdC4gQWxzbyBjbGltYXRlIGNoYW5nZSBpcyBsZWFkaW5nIHRvIGhpZ2hlciBkZW5zaXRpZXMgb2YgaGlnaC1lbGV2YXRpb24gZm9yZXN0LiAgCiAKIDQuIFRoZSBhbW91bnQgb2Ygc25vdyBhbmQgdGhlIHRpbWluZyBvZiBpdHMgbWVsdCBpcyBpbXBvcnRhbnQgZm9yIGEgbG90IG9mIHRoaW5ncyBsaWtlIGxpZmUsIGFsbW9uZHMsIGFuZCBzd2ltbWluZyBwb29scy4gIAogCgojIERhdGEgQXZhaWxhYmlsaXR5CgpXZSBhcmUgdXNpbmcgTmFzYSdzIDNtIGxpZGFyLWRlcml2ZWQgW0FTTyBkYXRhXShodHRwczovL25zaWRjLm9yZy9kYXRhL2Fzb18zbV9zZCkuIExpZGFyIGZpZ2h0cyBoYXZlIG9jY3VycmVkIHN0YXJ0aW5nIGluIDIwMTMsIGFsdGhvdWdoIHRoZXkgaW1wcm92ZWQgdGhlaXIgZXF1aXBtZW50IGluIDIwMTQsIHNvIGxldCdzIHN0YXJ0IHRoYXQgeWVhci4gVGhleSBoYXZlIGNvbmR1Y3RlZCBmbGlnaHRzIGFjcm9zcyBhIG51bWJlciBvZiBiYXNpbnMsIGJ1dCB3ZSdsbCBmb2N1cyBvbiB0aGUgTWVyY2VkIGFuZCBLaW5ncyBiYXNpbnMgc2luY2UgdGhlc2UgaGF2ZSBoYWQgdGhlIGNsb3Nlc3QgdG8gbmF0dXJhbCBmaXJlIHJlZ2ltZXMgYW5kIGxpa2VseSBtb3N0ICduYXR1cmFsJyBkaXN0cmlidXRpb25zIG9mIGdhcCB0eXBlcy4gSWYgd2Ugd2FudCBtb3JlIHRlbXBvcmFsIGFuZC9vciBzcGF0aWFsIHJlcGxpY2F0aW9uIHdlIGNvdWxkIGxvb2sgYXQgb3RoZXIgYmFzaW5zLCBidXQgd2Ugd291bGQgbmVlZCB0byB3cmFuZ2xlIGEgbG90IG9mIGRhdGEsIGFuZCBkYXRhIHJlcXVlc3RzIGhhdmVuJ3QgYmVlbiByYXBpZCBzbyBmYXIuICAgICAKCiFbU2llcnJhIE5ldmFkYSBiYXNpbnMgc3VydmV5ZWQgYnkgQVNPXShBU09fRmxpZ2h0X01hcC5ibXApICAKCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFfQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IEZBTFNFLCB3YXJuaW5nPUZBTFNFLCBtZXNzYWdlPUZBTFNFKQprbml0cjo6b3B0c19rbml0JHNldChyb290LmRpciA9ICIuLi8uLiIpCgpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShsdWJyaWRhdGUpCmxpYnJhcnkocGxvdGx5KQpsaWJyYXJ5KGNpdHIpCmxpYnJhcnkoa25pdHIpCmBgYAoKCkN1cnJlbnRseSB3ZSBoYXZlIGRhdGEgZm9yIDIwMTQgZnJvbSB0aGUgTWVyY2VkIGJhc2luLCBhbmQgYSByZXF1ZXN0IGluIGZvciBhbGwgb3RoZXIgZmxpZ2h0cyBpbiB0aGUgTWVyY2VkIGFuZCBLaW5ncyBiYXNpbnMuICAKCmBgYHtyfQpkIDwtIHJlYWQuY3N2KCJkYXRhL3JlZmVyZW5jZS9zbm93ZGVwdGgzbV9hdmFpbC5jc3YiKSAlPiUKICAgIG11dGF0ZShkYXRlID0gYXMuY2hhcmFjdGVyKGRhdGUpLAogICAgICAgICBkYXRlID0gYXMuRGF0ZShkYXRlLCAiJVklbSVkIiksCiAgICAgICAgIHllYXIgPSB5ZWFyKGRhdGUpLCAjeWVhcgogICAgICAgICBqZGF5ID0gYXMuaW50ZWdlcihmb3JtYXQoZGF0ZSwgJyVqJykpKQoKIyMgUGxvdCBmbGlnaHRzIGFjcm9zcyB5ZWFycyBhbmQgZGF5CmdncGxvdChkLCBhZXMoeCA9IGpkYXksIHkgPSB5ZWFyLCBjb2xvciA9IGJhc2luLCBzaGFwZSA9IGZhY3Rvcihwb2Nlc3MpKSkgKwogIGdlb21fcG9pbnQoc2l6ZSA9IDQpICsKICBzY2FsZV9zaGFwZV9tYW51YWwodmFsdWVzID0gYygyMSwgMTkpLAogICAgICAgICAgICAgICAgICAgICBsYWJlbHMgPSBjKCJSZXF1ZXN0ZWQiLCJPYnRhaW5lZCIpLAogICAgICAgICAgICAgICAgICAgICBuYW1lID0gIkRhdGEgU3RhdHVzIikgKwogIHNjYWxlX2NvbG9yX2JyZXdlcihwYWxldHRlID0gIkRhcmsyIiwKICAgICAgICAgICAgICAgICAgICAgbGFiZWxzID0gYygiS2luZ3MiLCAiTWVyY2VkIiksCiAgICAgICAgICAgICAgICAgICAgIG5hbWUgPSAiQmFzaW4iKSArCiAgeGxhYigiSnVsaWFuIERheSIpICsgeWxhYihOVUxMKQoKCmBgYAogCgojIFF1ZXN0aW9ucyAgCgpUaGVyZSBhcmUgKGF0IGxlYXN0KSB0d28gd2F5cyB0byBhcHByb2FjaCB0aGlzIGFuYWx5c2lzIHVzaW5nIDNtIEFTTyBkYXRhOiAgCgogIDEuIEF0IHRoZSBnYXAgbGV2ZWw6IEhvdyBkb2VzIGdhcCBzaXplIGFuZCBzaGFwZSBpbmZsdWVuY2UgbWVhbiBzbm93IGRlcHRoIGFuZCBzbm93IGRpc2FwcGVhcmFuY2UgZGF0ZT8gQ291bGQgYXBwbHkgdGhlIFNEQyBtZXRyaWMgaGVyZSBhcyB3ZWxsLiAgCiAKICAyLiBBdCB0aGUgcGl4ZWwgbGV2ZWw6IEhvdyBkb2VzIHNub3cgZGVwdGggYW5kIHByZXNlbmNlIHZhcnkgd2l0aCBkaXN0YW5jZSB0byBwYXRjaCBlZGdlPyBEbyBwaXhlbHMgaW4gc2hhZGVkIGVkZ2VzIGhvbGQgc25vdyBsb25nZXI/ICAKIApFaXRoZXIgd2F5IHdlJ2Qgd2FudCB0byBpbmNvcnBvcmF0ZSBlbGV2YXRpb24sIGFzcGVjdCwgZGF0ZSwgYW5kIHNub3cteWVhciB0byBhc2sgdGhlIGZvbGxvdyBpZiBwb3NzaWJsZTogIAogCiAzLiBEbyB0aGlzIHJlbGF0aW9uc2hpcCBkZXBlbmQgdXBvbiBjbGltYXRlIChib3RoIGludGVyLWFubnVhbCBhbmQgYWNyb3NzIGVsZXZhdGlvbik/ICAKIAoKIyBSZWxldmFudCBMaXRlcmF0dXJlICAKClRoZXJlIGFyZSBjbGVhcmx5IG1vcmUgdGhhbiBpcyBsaXN0ZWQgYmVsb3cuIElmIHRoZXJlIGFyZSBhbnkgaW1wb3J0YW50IG9taXNzaW9ucyBhdCB0aGlzIHBvaW50LCBwbGVhc2UgbGV0IG1lIGtub3cuICAKCgogLSBSZWxhdGVkIHRvIGZpcmUgZWZmZWN0cyBvbiBmb3Jlc3QgZ2Fwcy9wYXRjaGVzOiBbQHN0ZWVsQ2hhbmdpbmdMYW5kc2NhcGVXaWxkZmlyZTIwMTg7IEBzdGV2ZW5zQ2hhbmdpbmdTcGF0aWFsUGF0dGVybnMyMDE3OyBAY29sbGluc0FsdGVybmF0aXZlQ2hhcmFjdGVyaXphdGlvbkZvcmVzdDIwMTc7IEBib2lzcmFtZU1hbmFnZWRXaWxkZmlyZUVmZmVjdHMyMDE3XSAgCiAKIC0gTGlkYXIgbWVhc3VyZW1lbnRzIG9mIHNub3cgZGVwdGg6IFtAaGFycG9sZExpREFSZGVyaXZlZFNub3dwYWNrRGF0YTIwMTQ7IEBwYWludGVyQWlyYm9ybmVTbm93T2JzZXJ2YXRvcnkyMDE2XSAgCiAKIC0gU25vdyBhbmQgdHJlZSBjb3ZlcjogW0BoYXJwb2xkQ2hhbmdlc1Nub3dBY2N1bXVsYXRpb24yMDE0OyBAY3VycmllclNub3dEZXB0aFZhcmlhYmlsaXR5MjAxODsgQGx1bmRxdWlzdExvd2VyRm9yZXN0RGVuc2l0eTIwMTM7IEBzdGV2ZW5zU2NhbGVkZXBlbmRlbnRFZmZlY3RzUG9zdGZpcmUyMDE3XSAgCiAKIAoqKiogIAoKIyBCaWJsaW9ncmFwaHkgIAoKICAKCg==